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The  objectives  of  this  program  are  to:  j 

-  Evaluate  the  frequency  instabilities  of  quartz  crystal  resonators  ate  a 
function  of  temperature  and  Q- factors. 

-  Show  the  feaaability  of  low  temperature  quartz  oscillators 

-  Determine  the  dynamic  thermal  behavior  of  SAW  quartz  oscillators  and 

precise  crystal  orientation  and  wave  propagation  direction  with  minimized 
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>>  The  measurements  of  the  frequency  noisfev^f  different  quartz  resonators  at 
room  temperature  show  a  correlation  between T7ff  noise  levels  and  Q-factors, 
following  a  1/Q  law.  This  new  property  was  corff^r™*^  by  measuring  the  1/f 
noise  level  of  a  pair  of  resonators  at  300  K,  4.2?k\and  1  K. 

A  first  theoretical  study  enables  to  explain  the  1/Q*  law  from  the  phonon- 
phonon  interactions  due  to  the  crystal  anharmonicities. 

Under  fast  temperature  variations  frequency  shifts  are  observed  in  SAW 
delay  lines,  which  are  proportional  to  the  time  derivative  of  the  temperature. 
This  dynamic  thermal  behavior  was  evaluated  for  ST  cut  quartz  and  found  com¬ 
parable  to  that  one  of  AT-cut  quartz  BAW  resonators.  Therefore  such  a  thermal 
sensitivity  is  not  admissible  for  fast  warm-up  oscillator  applications  and 
must  be  lowered. 
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I  -  LOW  TEMPERATURE  QUARTZ  CRYSTAL  OSCILLATOR 


GENERAL  PRESENTATION 

The  frequency  stability  of  quartz  oscillators  is  affected  by 

-  the  noise  of  the  electronics  :  thermal  noise  and  1/F  noise 

-  the  peculiar  frequency  noise  of  the  resonator  :  frequency  random  walk 
noise  and  1/F  noise. 

When  the  noise  of  the  electronics  is  reduced  at  the  lowest  possible  level 
compatible  with  the  state  of  the  art  of  electronic  components,  the  limit  of  sta¬ 
bility  of  the  oscillator  corresponds  to  the  resonator  noise. 

-  the  frequency  random  walk  is  essentially  due  to  temperature  fluctua¬ 
tions  and  can  be  reduced  by  using  high  quality  temperature  controlled 

(2) 

ovens 

-  the  1/F  noise  level  indeed  defines  the  ultimate  limit  of  stability.  To 
the  1/F  spectrum  of  the  frequency  domain  corresponds  a  floor,  the  flicker  floor, 
in  the  time  domain. 

It  was  known  that  at  very  low  temperature  the  internal  losses  and  the 
acoustic  attenuation  of  quartz  material  are  reduced  and  therefore  the  Q-factor 
of  resonators  is  increased. 

The  present  work  shows  that 

-  the  i/F  noise  level  at  room  temperature  is  correlated  with  the  Q-factor, 

4 

following  a  1/Q  -law. 

-  at  low  temperatures  it  is  possible  to  minimize,  by  at  least  one  order 
of  magnitude,  the  1/F  noise  level.  Improvement  of  the  frequency  stability 
follows . 


1-1 


I 


The  study  of  the  noise  levels  as  a  function  of  temperature  required 

-  cryogenic  equipment  operating  at  1°K 

-  test  set  for  crystal  adapted  to  very  high  Q's  and  low  motional  resis¬ 
tances 

-  frequency  noise  measurement  system  for  pairs  of  resonators  at  300°K  am; 
1  °K . 


1/F  noise  is  an  universal  noise,  which  can  be  observed  in  almost  all  phy¬ 
sical  phenomena  and  systems.  Its  origin  and  mechanisms  are  not  known.  A  theore- 

4 

tical  study  has  been  undertaken,  and  the  first  task  was  to  explain  the  1/0  de¬ 
pendence  law,  which  was  experimentally  demonstrated. 


CRYOGENIC  EQUIPMENT 

A  liquid  helium  dewar  with  double  isolation  was  used.  Its  capacity  is  5 
and  gives  an  autonomy  of  one  or  two  days,  depending  on  the  number  of  leads  and 
cables  sunk  in  it.  This  small  size  was  choosen  in  order  not  to  expend  to  much 
liquid  helium  when  the  different  quartz  samples  successively  are  cooled  and 
heated.  In  the  second  phase  of  the  project  a  30  liters  dewar  will  be  used  to 
maintain  Xtal  oscillators  below  4K  during  periods  of  time  of  the  order  of  seve¬ 
ral  weeks.  Vacuum  pumps  enable  to  reach  temperatures  between  4.2K  and  IK.  Tempe 
rature  is  measured  by  means  of  carbon  and  germanium  resistors,  and  also  by  mean; 
of  the  bath  pressure  when  below  4.2K. 

For  guide  measurements,  the  resonators  are  directly  sunk  in  the  helium 
bath.  If  more  stability,  mainly  versus  temperature,  is  required  they  are  enclo 
sed  in  a  copper  can. 

This  cryogenic  equipement  was  ready  to  use  in  February  1980. 
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CRYSTAL  PARAMETERS  AT  l,OW  TEMPERATURE 


The  measurements  were  performed  by  using  the  *  transmission  method. 


Measurement  system 

Measurements  at  low  temperature  present  some  difficulties  because  of  the 
length  of  the  cables  between  the  resonator  at  the  bottom  ol'  the  dewar  and  the 
source  and  detector  which  are  at  room  temperature.  Therefore  it  is  necessary 
to  place  the  *  network  with  the  resonator  in  the  bath. 

The  test  set  is  represented  on  fig.  1. 
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The  n  network  must  be  adapted  to  50  Cl  ,  the  characteristic  impedance  of 
the  cobles.  But.  the  vnlues  of  the  parallel  resistors  must  be  chooses  as  ,i  (uiw- 
tion  of  the  motional  resistance  R^  of  the  resonator.  This  last  one  can  become 
close  to  1  Cl  at  very  low  temperature.  Therefore  two  »  networks  were  used.  Theii 
values  are  given  on  fig.  2a  and  2b. 


The  n  network  head  must  be  constructed  with  care.  To  avoid  spurious  con¬ 
tact  resistances  at  temperatures  below  4K  all  the  electrical  contacts  were  gold 
platted.  The  n  head  is  represented  on  fig.  3. 


pins  to 
connect 
Xtal 


Frequency,  motional  resistance  and  Q-factor  were  measured  between  IK  and 
20K.  This  temperature  range  corresponds  to  the  diminution  of  the  acoustic  atte¬ 
nuation,  the  phonon-phonon  interaction  peak  being  located  at  20K, 

The  highest  Q-factor  was  obtained  with  G4  and  corresponds  to  Q  =  150.10^ 
at  T  =  1 .5K. 
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and  5. 


The  Q  and  R  values  of  the  resonators  provided  by  RADC  are  shown  on  fig.  4 

The  corresponding  frequency  shifts  are  given  on  fig.  6.  They  are  of  the 
order  of  5  to  6  KHz.  the  temperature  coefficient  of  the  frequency  was  also  meac 
red.  It  decreases  with  temperature,  as  expected  (fig.  7). 

At  room  temperature,  the  O-factor  values  of  these  resonators  ar<>  clos'd, 
hut  at  low  temperature  a  large  dispersion  is  observed.  The  inverse  O-factor  can 
be  written 

1/0  =  1/0.  1/0  +  1/0 
b  s  e 

This  relation  points  out  the  different  loss  contributions  due  to  the  crystal 

bulk  (1/0,  ),  the  surface  (1/0  )  and  the  environment  (1/0  )  constituted  of  the 
b  s  e 

electrodes  and  the  mounting.  1/Q^  is  representative  of  the  crystal  properties 
and  depends  on  the  phonon  interaction  processes.  Its  contribution  is  predominan' 
at  20°K  and  certainly  also  at  room  temperature.  1/Q^  decreasing  with  temperatur' 
as  T4  it  becomes  negligible  at  low  temperature.  The  influence  of  1/Q^  and  1/Q^ 
appears  and  leads  to  a  threshold.  These  two  last  contributions  are  mainly  repre 
sentative  of  the  resonator  technology. 


The  1/f  frequency  noise  of  resonators  cannot  be  measured  in  an  oscillatin/ 
loop,  because  the  phase  fluctuations  of  the  loop  due  to  the  electronics  are  gent 
rally  higher.  The  method  which  is  used  consists  in  driving  in  a  it  transmission 
network  two  identical  resonators  by  means  of  a  frequency  source  with  high  spec¬ 
tral  purity  (3) (4).  j<jentical  resonators  means  resonators  with  identical  frequer 
cies  and  0-factors,  or  made  identical  with  tuning  capacitors  and  resistors.  The¬ 
refore  a  preselection  of  pairs  of  resonators  at  the  same  frequency  is  necessary. 
The  measurement  system  is  shown  on  fig.  8. 
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When  the  two  arms  of  the  phase  bridge  are  correctly  balanced  SO  dB  to  60 
iJB  of  rejection  of  the  source  frequency  fluctuations  are  obtained. 

The  resolution  depends  on  the  resonators  0-factors  following  the  relation 
between  the  phase  and  frequency  fluctuations  d+  =  -  2Q  df/f  ,  for  Fourier  fre¬ 
quencies  within  the  resonator's  bandwidth. 

On  fig.  9  is  shown  a  typical  frequency  noise  spectrum.  Three  different 
noises  appear  :  frequency  random  walk  (1/F1)  for  the  lowest  Fourier  frequencies 
flicker  noise  (1/F)  for  the  intermediate  ones  and  white  noise  (F°)  for  the 

highest  ones.  1/F1  noise  is  related  to  the  crystal  temperature  fluctuation  AT 

^  ( 3 ) 

through  its  dynamic  thermal  coefficient  a  following  the  law  Af/f  =  a  AT/At 

This  noise  can  be  minimized  by  using  thermostats  with  high  accuracy.  For  AT  cut 

\  _13 

a  is  10  s/°C.  This  means  that  the  achievement  of  a  10  frequency  stability 

—8 

needs  a  temperature  control  as  good  as  10  °C/s. 

The  white  noise  which  is  observed  has  nothing  to  do  with  the  resonators, 
but  is  due  to  the  measurement  system  itself. 


The  level  of  1/F  noise  was  measured  for  a  large  number  of  pairs  of  resonn 
tors  .it  different  frequencies,  and  with  a  large  dispersion  of  the  0-factor  va¬ 
lues.  On  fig.  10  are  plotted  the  different  noise  levels  measured  on  all  those 
resonators  and  given  at  1  Hz  from  the  carrier. 


From  these  experimental  data  the  dependence  between  1/F  noise  and  Q-facto- 
is  determined  by  using  a  regression  method  (linear  regression  with  the  loganth 
mic  scales).  This  leads  to  the  general  relation 


The  values  obtained  for  n  and  1  depend  on  the  selection  of  the  experimen¬ 
tal  points  among  the  18  available  values.  After  mathematical  treatment  it  can 
be  shown  that  3.8  4  n  ■?  4.3.  Therefore  it  is  proposed  to  consider  that  the  most 
probable  value  will  be  n  =  4 .  A  small  variation  in  the  n  value  induces  a  drasti 
change  of  the  constant  X,  and  for  the  previous  range  of  n  the  corresponding  va¬ 
lues  of  x  are  0.1  <  X  <  60.  A  reasonable  fitting  will  correspond  to  x  close 

to  unity  when  n  is  made  equal  to  the  integer  4. 


1-12 


2 )  Low  temperatures 


At  room  temperature  there  is  no  way  for  modifying  the  Q-factor  of  one  re¬ 
sonator.  An  additional  resistor  will  change  the  loaded  0,  but  not  the  unloaded 
one  which  is  relevant  for  the  noise  measurements.  The  only  possibility  was  to 
use  different  resonators  with  different  unloaded  Q-factors. 

But  the  acoustic  attenuation  is  strongly  dependent  of  temperature  at  low 
t  empe-'ature;:  and  ran  be  either  increased  near  the  phonon-phonon  interaction  pea 
and  near  the  impurity  relaxation  peaks,  or  decreased  for  temperatures  below  10K 
as  it  has  been  shown. 

The  previous  noise  measurement  system,  which  requires  a  pair  of  resonat'- 
at  the  same  frequency,  is  not  aaapted  for  low  temperature  measurements.  The 
frequency  variations  from  300K  to  4K  are  of  the  order  of  a  few  KHz.  But  because 
of  small  differences  in  the  crystal  orientation,  this  frequency  variation  vers- 
sus  temperature  is  not  reproductible  enough.  Even  when  identical  at  room  tempe¬ 
rature,  the  frequencies  of  a  pair  of  resonator  present  differences  of  several 
100  Hz  which  cannot  be  cancelled  out  by  using  tuning  external  capacitors. 

Therefore  a  new  measurement  system  for  low  temperatures  was  studied.  Its 
schematic  diagram  is  represented  on  fig.  11. 


F*f  ISO 


(•'!«.  11 
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Each  of  the  two  resonators  is  excited  by  one  of  the  side  bands  of  an  ampl 
tude-modulated  signal  (without  carrier).  At  the  output  of  the  first  mixer  the 
intermediate  frequency  is  affected  by  the  total  frequency  noise  of  the  resona¬ 
tors.  The  influence  of  the  driving  source  is  cancelled  out  by  balancing  the  cha 
racteristics  of  the  resonators.  It  is  necessary  to  use  a  good  frequency  source 
for  the  modulating  signal,  even  if  its  frequency  spectrum  in  purified  by  divi¬ 
ding  by  4000 . 

The  first  mixer  is  followed  by  a  low  noise  amplifier,  with  20  dB  of  gain. 

A  second  mixing  with  a  2f  reference  frequency  delivered  from  a  digital 
phase  shifter  coupled  to  the  divider  gives  a  noise  voltage  proportional  to  the 
resonator's  frequency  fluctuations.  A  filter-amplifier  cancells  out  the  4f  fre¬ 
quency  component.  A  FFT  spectrum  analyzer  gives  directly  the  S^(F)  spectrum  of 
the  two  resonators. 

One  other  advantage  of  such  a  measurement  system,  when  compared  to  the 
previous  one,  is  to  enable  easier  and  faster  calibration,  just  changing  by  a 
small  amount  the  modulation  frequency  and  measuring  the  DC  voltage  shift  at  the 
mixer  output. 

The  1/F  noise  of  a  pair  of  identical  AT-cut  quartz  resonators  (manufactu¬ 
red  by  CEPE)  were  measured  at  3 00K,  A.2K  and  at  IK.  The  frequency  difference 
between  the  two  resonators  waB  1  KHz.  The  noise  levels  for  one  resonator  are 
given  on  table  I. 


Temperature 

300  °K 

4 , 2°K 

1°K 

Q 

2.4. 106 

5.106 

9.106 

S  (1  Hz) 
y 

1.6.10-24 

1.4.10-25 

L  - 

3.5.10-26 

Table  I 


These  results  confirm  the  dependence  of  1/F  noise  on  Q-factor  and  in  this 

4 

case  with  the  same  resonator.  The  data  does  not  follow  strictly  the  1/0  law. 

A  comment  on  this  point  will  be  given  in  the  next  section. 
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theoretical  approach 


It  is  known  that  the  acoustic  attenuation  results  from  the  interaction  of 
the  sound  wave  with  the  thermal  phonons  of  the  crystal.  This  is  a  fundamental 
phenomenon,  which  would  appears  even  in  a  perfect  crystal,  i.e.  without  defects, 
because  of  the  anharmonicities  of  the  lattice,  due  to  the  nonlinear  nature  of 
the  interatomic  bonds.  Also  because  of  these  nonlinearities  are  modifications 
of  the  phonon  frequencies  and  of  the  sound  wave  velocity. 

Calculations  of  attenuation  were  made  by  using  two  different  theoretical 

( 5) 

approaches.  In  the  Landau  and  Rumer  theory  (1937)  the  sound  wave  is  conside¬ 
red  as  a  beam  of  phonons  which  are  scattered  by  collisions  with  the  thermal  pho¬ 
nons.  This  theory  does  not  take  into  account  the  interactions  between  thermal 
phonons  and  is  valid  only  in  the  low  temperature  range,  when  the  phonon  life 
time  and  mean  free  path  are  large  and  therefore  energy  and  momentum  are  well 
defined.  At  higher  temperatures  the  life  time  becomes  shorter.  The  uncertainties 
of  energy  and  momentum  increase  and  the  selection  rules  arising  from  energy  and 

momentum  conservation  laws  break  down.  In  the  second  theory  due  to  Akhieser 

( 6 ) 

(1939)  the  sound  wave  is  represented  macroscopically  (and  not  microscopi¬ 
cally  as  in  the  Landau-Rumer  theory).  The  thermal  phonon  system  is  disturbed 
from  equilibrium  by  the  sound  wave  strain  but  tends  to  return  to  equilibrium 
because  of  the  collisions  between  thermal  phonons.  In  this  case  the  phonons  are 
localized  and  their  mean  free  path  has  to  be  small  when  compared  to  the  sound 
wave  length. 


1 )  Harmonic  approximation 

The  main  results  of  the  lattice  dynamics  in  the  harmonic  approximation 
will  be  briefly  presented,  (for  more  details  cf.  ref.  7). 


The  potential  energy  $2  is  assumed  to  be  a  quadratic  function  of  the  ato¬ 
mic  displacements  uion<  where  n  refers  to  nth  Wigner-Seitz  cell,  a  to  the  ath 
atom  in  the  unit  cell  and  i  distinguishes  the  three  components  of  the  displace¬ 
ment  vector 


l 

‘  ian 


i'a'n' 

6 

vian 


uian  ul'o'n' 


(1) 


i  *  o  *  n  ' 

are  the  second-order  interactomic  coupling  parameters. 

ian 
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The  kinetic  energy  of  the  crystal  is 


V  =  (*>  [  (2) 

ian  d 

P.  is  the  ith  component  of  the  momentum  of  the  ath  atom  in  cell  n,  and  M  is 
ion  ° 

the  mass  of  this  atom. 

The  Hamiltonian  of  the  crystal  is  then 

h2  =  +  *  (3> 


The  displacement  u.  and  momentum  P.  can  be  written  as  a  sum  of  waves 
ian  ian 

with  wave  vector  k,  polarization  j,  travelling  in  the  x(l)  direction. 

Vi  2nik.x(n) 

Uion  =  J.  [2NM>TkJ)]  eia(kJ)  I  e  (4) 

.  „  ..  %  2»  ik.x(n) 

P.an  =  -i  [  [  hMa^^i]  eia  { k  j )  [  a(kj)+a+(-kj)]  e  (5) 

kj 


is  the  polarization  vector,  N  the  total  number  of  unit  cells.  a+(kj)  and 
o(kj)  are  the  creation  and  annihilation  operators. 


It  can  be  shown  that  the  Hamiltonian  takes  the  simpler  form 

H  =  l  ti  u  ( k  j )  [  a+  ( k  j  )  a  ( k  j  )  +  %  ]  ( 6 ) 

kj 

The  sum  over  k  is  extended  over  the  first  Brillouin  zone. 


The  eigenvalues  of  the  Hamiltonian  H2  gives  the  energy  of  the  crystal 

E  =  [I  N(kJ)  ♦  '4]  n  u(kj)  (7) 

kj 

which  can  be  interpreted  as  the  sum  of  the  partial  energies 

E(kj  )={  N(kj  )+/$]tiw(kJ )  of  each  mode  k  and  polarization  j.  N(kJ)  is  an  integer.  L 
teem  of  phonons  it  indicates  the  number  of  phonons 'with  energy  "h*)  (kj )  in  each 
partial  wave. 
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2 )  Anharmonicity 


tn  the  harmonic  approximation  the  various  lattice  waves,  or  phonons,  are 
independant  and  there  is  no  transfert  of  energy  between  them.  Therefore  a  sound 
wave,  which  corresponds  to  a  given  lattice  wave  with  a  large  number  of  phonons, 
will  not  be  attenuated. 

a )  Landau-Rumer  theory 


Anharmonici ties  can  be  introduced  by  means  of  third  (and/or  fourth)  order 
interatomic  coupling  parameters.  In  quantum  mechanical  model  the  anharmonic  con 
tribution  to  the  Hamiltonian  takes  the  form 


H 


3 


3/2 

— H— .  y  _  1  (-kj,k'j',k"j») 

3/p  y  L  - - — - - - 

2  '  6N  7  kk'k"  A>(kj)  wfk"^') 

j  j  '  j" 


a+(kj )  a(k' j 1 )  a(k" j" ) 


+  all  other  possible  combinations  of  the  a  and  a+  operators. 


(a) 


These  three  phonon  interactions  must  verify  the  energy  and  momentum  consei 
vation  laws. 


The  sound  wave  can  be  considered  as  a  lattice  wave  of  mode  KJ  in  which 
phonons  are  added  so  that  their  number  N(kj)  ia  in  excess  of  the  thermal  cquili 
brium  number  n(kj).  The  rate  at  which  phonons  are  scattered  out  of  thi3  mode 
enables  to  define  a  lifetime  t(KJ)  given  by 


1  _  _ 1 _  dN(KJ) 

1 ( KJ )  "  “  N (KJ )  -  n(KJ )  dt 


(9) 


If  s(KJ)  is  the  phase  velocity  of  the  wave  the  attenuation  a  is  a  =  \/z  is. 


The  rate  of  change  dN(KJ)/dt  can  be  calculated  by  using  the  Golden  Rule 
of  the  perturbation  theory. 


The  anharmonic  coupling  also  perturbs  the  energy  levels  of  the  sound  wave. 
This  perturbation  is  calculated  by  using  the  standard  second-order-time-lndcpen- 
deiTt  perturbation  theory. 

The  validity  of  this  method  can  be  evaluated  by  considering  thnt  the  the’- 

mal  phonons  have  a  finite  life  time.  Let  t  be  lte  average.  In  account  of  the 

Uncertainty  Principle  this  finite  lifetime  makes  the  phonon  energy  uncertain 
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by  an  amount  of  the  order  of  h/t.  Therefore  the  method  will  not  be  valid  as  soor 
as  the  energy  uncertainty  is  comparable  to  the  sound  wave  phonon  energy  lift.  Thu 
the  Landau-Rumer  method  is  valid  when  ffx  «  lift,  or  {It  >>  1.  This  condition  is 
satisfied  for  sound  waves  at  high  frequencies  or  at  low  temperatures. 


b)  Akhieser  Theory 

At  higher  temperatures  the  thermal  phonon  lifetime  decreases.  As  a  conse¬ 
quence  their  mean  free  path  becomes  smaller  than  the  sound  wave  length  X.  The 
phonons  can  be  considered  as  localized  wave  packets.  The  thermal  phonon  avera¬ 
ged  frequency  at  temperature  T  is  k  T/1S  (k_  :  Boltzmann  constant).  If  l  is  the 

D  B 

wave  packet  length,  the  condition  1  <<X  is  equivalent  to  "tin  <<  knT. 

D 

In  this  case  the  method  proposed  by  Akhieser  consists  in  using  a  macros¬ 
copic  model  for  the  sound  wave  and  the  microscopic  one  for  the  thermal  phonons. 
The  anharmoniclties  are  introduced  by  means  of  a  strain  dependance  of  the  phonoi 
frequencies.  A  detailed  description  of  the  method  can  be  found  in  ref  (8). 


The  total  energy  of  the  crystal  of  volume  V  is  the  sum  of  the  energy  of 
the  thermal  phonons 

E  =  l  (N(kj)  +  !0tl  w  (kj )  (10) 

kj 

and  of  the  mechanical  potential  energy  due  to  the  strains  induced  by  the  sound 
wave 

E  .=C„Ti-+JiC„,n„n.  +  ...  (11) 

sound  aB  aB  aByfi  aB  oy 

C  „  and  C  „  ,  have  the  dimension  of  elastic  constants. 

OB  GtpYti 

The  derivative  of  the  total  energy  with  respect  to  strain  gives  the  stress 


T  =  ( a /a 

aB 


n  ) 

aB 


th 


+  E  .) 
sound 


(12 


"  When  in  equilibrium  the  number  of  phonons  in* mode  kj  is  given  by  the  Bose- 
Einstein  distribution 


n(kj)  =»  (exp  [Tj  w(kj)  /  kgT]  -  l)"1 


(13) 
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and  this  equilibrium  corresponds  also  to  T  =  n  „  =  0. 

aB  aB 


Inis  enables  to  determine  completely  the  stress 


TaB  =  CaBy*  \«  +  J  I  tN<kJ>  " 


kj 


kj  aB  y« 


an 


aB 


(Id) 


The  mechanical  equilibrium  equation 


B 

— r 

a  t 


a  t 


aB 


ax 


and  the  strain-displacement  relation 


au 


a  u 


n  -  ♦  — — 

aO  ?  '  3x  3 x 
B  a 


(15) 


(16) 


enables  to  calculate  plane  waves  solutions  with  wave  vector  K  and  frequency  .fl¬ 
it  2*K.x-nt) 


u„  =  u  e  (KJ)  e 

B  o  B 

(17) 

i(2»K.x-nt) 

N (k j )  =  n(kj)  +  iN(kj)  e 

(IS) 

In  the  absence  of  phonons  the  solution  has  the  simpler  form 

,  (  i(2*K.x-Sl  t) 

u  =  u  e  (KJ)  e  0 

B  o  6 


(19) 


The  correction  due  to  the  phonons  can  be  calculated  by  using  a  perturba¬ 
tion  method.  This  leads  to 


”  ■  "o  *  |5r  V«  J  t kj j 

o  kj  aB  y6 

_  igL  e  (KJ)  K  i  mm 

pVf!o  °  8  kj  3naS  Uo 


(20) 


The  attenuation  a  is  given  by 

n 


a  =  -  Im 


3  (KJ) 
o 


(21) 
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and  the  new  velocity  of  the  sound  wave 


These  quantities  are  completely  determined  when  AN(kj)  is  known.  This  cor 
responds  to  the  calculation  of  the  response  of  the  phonon  system  to  the  sound 
wave,  which  is  given  by  the  Boltzmann  equation 


3N(kj)  1  3o.(k,j)  3nae  3N(  k  j )  3N(k,j)  ,  3N(kj)  , 

3 1  "  2*  3n  „  3x  3k  _  V  3x  +  1  3t  J  coll 

aB  Y  Y  a 


Remark  :  The  term 
-  3t  coll 

the  collisions  between  phonons. 


is  most  important  because  it  takes  into  account 
Without  collisions  and  with  tlr  >>  1  and 


n  <<  k^T/fl  the  present  theory  goes  back  to  the  Landau-Rumer  theory. 


Three  different  types  of  collisions  can  be  distinguished  : 

-  Normal  processes  (N),  in  which  quasimomentum  and  energy  of  thermal  pho¬ 


nons  are  conserved  in  three  phonon  collisions 

hk  =  hk'  +  hk"  (24) 

flu)  =  ttu'  +  flui"  ( 25 ) 

Total  energy  and  total  quasimomentum  also  are  conserved.  Therefore 

l  ■»> 

l  ( 3N(kj )/at  1  hk  =  0  (27) 

kj  C0X1 

The  result  of  N  processes  is  a  distribution  law  which  has  the  form 

[ftui(kj)  -  -nvk]/k  T 

NN(kj)  =  (e  0  -U-1  (28) 


where  7  and  T  are  such  that  the  conservation  equations  are  satisfied. 

-  Umklapp  processes  ( U ) 

The  quas i momentum  is  not  conserved 

hk  =  hk'  +■  hk"  +  hg  (29) 

•flu>  =  flui*  +  flw" 
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(30) 


In  the  case  of  U  processes,  of  N  processes  with  some  E  and  U  processes 
(necessary  to  give  a  redistribution  of  the  quasimomentum),  or  of  E  processes 
with  some  N  and  U  processes  (also  to  bring  the  phonon  distribution  to  equili¬ 
brium)  attenuation  and  velocity  take  the  general  form 


(35) 

(36) 


where  y  is  an  effective  Grueneisen  constant  depending  on  the  nature  of  the  pro¬ 
cess  involved,  C  the  specific  heat  ;  and  t  the  mean  relaxation  time  of  the  ther¬ 
mal  phonons . 
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Introducing  the  low  frequency  young  modulus  Eq  =  os^,  the  quantity  AE=CTtj 

and.  the  relation  between  attenuation  and  0-factor  a  =  fl  /(2Qs  ),  and  between 

o  o 

velocity  and  frequency  (s-3  )/s  =  (tl-fi  )/Si  one  obtains 

oo  oo 


1 

0 


AE 

E 

o 


n 


o 


T 


1+n’T1 

o 


(37) 


(9) 

In  the  case  of  a  resonator  fluctuation  4r  in  the  relaxation  time  t  will  pro¬ 
duce  a  fractional  frequency  fluctuation  y  *  6(l/fto  the  resonator  frequency 


y 


AE 

E 

o 


n'-t 

o 

(l+fl't 

o 


s 


*) 


6  T 
T 


(39) 


In  terms  of  power  spectral  densities  S^(F)  and 
(37)  in  (39)  one  obtains 


( F)  and  using  the  relation 


.4  'aE' 


5«t/t 


(40) 


4 

This  relation  shows  the  1/0  law  which  was  observed  experimentally,  bet¬ 
ween  the  noise  levels  and  the  Q-factors. 

For  higher  values  of  fJQT  the  Landau-Rumer  method  must  be  used.  At  the  pre¬ 
sent  state  of  this  study  the  frequency  change  at  low  temperature  was  not  calcu¬ 
lated  yet.  Only  attenuation  was  evaluated. 

Near  room  temperature  attenuation  is  proportional  to  T  as  shown  by  rela- 

4 

tion  (35).  At  low  temperature  it  becomes  proportional  to  T  .  This  corresponds 
to  what  is  observed  experimentally,  even  if  the  attenuation  behavior  is  someti- 

4  6 

mes  closer  to  T  or  T  .  The  difference  can  be  explained  by  the  temperature  depen 
dance  of  t ,  which  does  not  appears  explicitely  in  the  relations. 

However  temperature  also  appears  in  relation  (40)  in  both  0  and  E.  When 

T  decreases,  0  increases  and  therefore  there  is  some  compensation  in  the  product 
4  2 

0  E  .  This  can  explain  why  the  three  experimental  points  of  table  I  correspon¬ 
ding  to  1/F  noise  levels  at  300°K,  4.2K  and  IK  does  not  strictly  follow  the  1/0 
law. 
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CONCLUSION 


The  experimental  results  obtained  during  this  study  shown  for  the  first 
time  a  correlation  between  1/F  noise  and  O-factor  in  quartz  crystal  resonators. 

This  dependence  was  confirmed  on  the  one  hand  by  low  temperature  measure¬ 
ments  and  on  the  other  hand  by  a  theoretical  model. 

This  theory  based  on  phonon  interactions  is  a  first  approach  to  a  more 
general  theory  of  flicker  noise.  For  this  purpose  it  will  be  necessary  to  pre¬ 
cise  the  fluctuation  spectrum  of  the  phonon  relaxation  time,  and  calculate  the 
frequency  shift  at  low  temperature  by  using  the  Landau-Rumer  theory,  in  order 
to  introduce  completely  the  temperature  dependence  in  the  model. 

For  a  practical  point  of  view,  the  large  increase  of  Q-factor  and  the  lo¬ 
wer  noise  levels  open  the  possibility  of  achieving  cryogenic  quartz  crystal  os¬ 
cillators  with  better  stability. 
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II  -  FAST  WARM-UP  SAW  OSCILLATOR 


GENERAL  PRESENTATION 


The  characteristics  of  surface  acoustic  waves  are  modified  by  temperature 
variations.  Those  external  effects  are  at  the  origin  of  the  main  causes  of  fre¬ 
quency  instabilities  in  HAW  oscillators. 

When  the  temperature  fluctuations  are  slow  enough,  the  thermal  distribu¬ 
tion  in  the  volume  of  the  crystal  is  uniform.  The  usual  model  using  a  polynomia 
expansion  up  to  the  third  order  versus  temperature  is  satisfactory  to  describe 
the  static  thermal  behavior.  For  faster  temperature  variations  a  spatial  distri 
bution  of  temperature  within  the  crystal  arises.  Consequently,  thermal  effects 
involve  temperature  gradients  which  induce  thermal  stresses  and  strains.  This 
is  the  "Dynamic  thermal  behavior" ^  ^ ^ ^ .  These  thermal  stresses  and  strains  arc 
coupled  with  the  surface  wave  by  the  intrinsic  and  induced  nonlinearities  of 
the  substrate.  Velocity  and  frequency  shift  follow. 


The  present  work  shows  that 


-  the  fast  temperature  variations  induce  frequency  changes  which  are  pro¬ 
portional  to  the  time  derivative  of  the  internal  temperature. 


-  the  measured  ST-cut  dynamic  temperature  coefficient  a  is  equal  to 

-5  ST 

2.10  s/°K  and  is  comparable  to  the  coefficient  of  the  AT-cut  quartz  crystal 

—5 

bulk  wave  resonator  (1.3.10  s/°K). 


i  ' 

-  the  measured  dynamic  TC's  of  doubly  rotated  cuts  are  larger  in  SAW  than 


in  BAW  oscillators. 


The  theoretical  study  was  composed  of 

-  calculation  of  the  temperature  distributions  and  corresponding  thermal 
stresses  and  strains  by  taking  into  account  the  mechanical  mounting  of  the  SAW 
device . 

<V 

-  calculation  of  the  dynamic  T.C.  a  by  using  a  perturbation  method  applle. 
to  the  general  model  of  an  acoustic  wave  propagating  in  a  strained  medium. 
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The  measurements  of  the  dynamic  thermal  behavicr  required  a  temperature 
controlled  oven  operated  by  means  of  a  computer.  Numerical  programs  enable  to 
perform  given  temperature  cyclings,  from  the  corresponding  frequency  responses 

M 

the  a-coeff ic ient  are  deduced. 


TEMPERATURE  DISTRIBUTION 


A  unidimensional  model  has  been  used.  It  consists  in  a  crystal  plate  of 
2h-thickness  and  of  infinite  lateral  sizes  along  a^  and  (fig.  1). 


*(t) 


Fig.  1  :  Temperature  distribution 


A  more  accurate  model  should  take  into  account  the  influence  of  the  ford 
at  the  fixation  points  under  the  crystal  thermal  expansion.  Only  the  bottom  of 
the  plate  was  submitted  to  temperature  variations  #(t). 

Let  Gfa^.t)  be  the  temperature  inside  the  crystal.  The  boundary  condition 

are 

0(-2h,t)  =  *(t)  and  6(0, t)  =  0  (1) 

The  heat  transfer  with  the  surrounding  medium  on  the  propagation  surface 
(»2  =  0)  ia  considered  to  be  very  fast. 

Consequently,  temperature  diatribution  within  the  crystal  follows  the  gen 

(3) 

ral  relation  given  by  H.  Carslaw 
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2  2  , 
n  ir  kX 


->  (-1)“  $(X)dX 


where  K  is  the  thermal  diffusivity  constant  of  the  crystal  in  the  direction. 

Considering  that  «(t)  is  a  slowly  varying  function  of  time,  the  integra¬ 
tion  of  I(t)  has  been  performed  by  parts,  and  the  time  derivatives  higher  than 
the  first  order  neglected. 


Then,  the  temperature  distribution  is  given  by 


a  1  /  a  \ 

O(a  t)  -  -  -1  *(t)  +  —  (  -  —  +  a  h  )  &(t) 

2h  3<  \  4h  1  ) 


$(t)  is  the  first  time  derivative  of  ♦ ( t ) 


*(t)  =  ♦(t),t 


The  internal  temperature  0(a_,t)  having  a  spatial  distribution,  thermal 

*  (4) 

stresses  and  strains  appear.  By  following  Holland  ,  the  adiabatic  deformatio. 
n^j  which  takes  place,  is  given  by  the  following  relations  written  in  the  matri 
cial  reference  system. 

(^]=  AN  a2  +  BN  (6) 


with  N  s  1,  3,  5 


[nz)  =  e(a2,t)  („’]  -  [C ZZ)-X  [CZN)[nN] 


with  Z  =  2,  4,  6 


The  different  coefficients  are  given  by  the  relations 


*»  ■  li?  '°h  *  e» 


2-3 
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B - —  (6  ~  +  ~) 

N  (2h)’  3  2 


KJ  -  *  ^ZZ^  ^zJ  ^ 


0 

o  =  /  9(a  ,t)  da 

-2h 


0 

6=1  0(a2,t)  a2  da2 

-2h 


(10) 

(11) 

(12) 


where  and  are  the  thermal  dilatation  constants  and  [C]  are  the  elastic 
constants . 


By  using  the  stress-strain  relations 
mations  and  the  temperature  distribution, 
reference  system) 


,  tensions  t. . 

ij 

are  calculated 


induced  by  the  defor- 
(also  in  the  matricial 


'ltN  '  ~l^CNN  '  “  ‘CNZ^CZZ 
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tczNl> 
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0(a2,t)  ) 


(13) 


and  [t^l  =  °  (14) 

in  account  of  the  boundary  conditions  corresponding  to  a  plate  whic 
is  free  to  expand. 


FREQUENCY  SHIFTS 


Considering  the  temperature  variations  as  a  bias  applied  to  the  crystal, 
propagation  properties  are  modified  and  in  a  natqral  state  coordinate  system 
the  nonlinear  propagation  equation  is  written 


po  ui,tt 


(AikJm  Uj,m),k 


(15) 


where  p  i s 

o 

vibrations. 


the  specific  mass  and  u 


J 


the  displacement  due  to  the  high  frequency 


The  boundary  conditions  corresponding  to  a  stress  free  surface  are 


A.,  . 

ikjm 


u  . 

Jim 


0 


=  0  for 


(16) 


where  A.,  .  are  the  modified  elastic  constants,  which  can  be  written  in  the  for 
ikjm 


A  =  C . .  .  +■)!,,. 

ikjm  ikjm  ikjm 


(17) 


where  the  H..  .  tensor  is  considered  as  a  small  term  with  respect  to  the  se 
ikjm 

cond  order  elastic  constant  C.,  .  .  It  can  be  related  to  the  temperature  and  the 

ikjm 

thermal  stresses  and  strains  in  the  following  way 


H  .  =  6  t,  +  U.  C,  .  +  U  C,  .  +  h  C.  +  a„  C . .  .  0(  a  ,  t )  ( 1R 

ikjm  ij  km  i,p  kpmj  j,q  kimq  uv  kmpquv  Cikjm  lkjm  2 


where  t,  ,  n  and  V,  are  respectively  the  thermodynamic  tensions,  the  defor- 
km  uv  i,p 

mations  and  the  displacement  gradients.  Here,  a ^  is  the  first  order  tempera 

ikjm 

ture  coefficient  of  elastic  constants  for  the  considered-cut. 

In  the  case  of  SAW  oscillators,  the  nonlinear  coupling  of  the  crystal  witl 

the  high  frequency  wave  depends  on  its  depth  penetration.  Then  perturbation 

terms  H~  are  function  of  the  step  variable  a„  and  can  be  written  following 
ikjm  2 

the  relation  : 


- -  r  (n)  n 

ikjm  L  ikjm  2 


(19) 


(5) 

By  using  a  perturbation  method  ,  relative  frequency  shifts  are  calcu¬ 


lated 


l 
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A,  v<’>  a;  •,;«>•  ■£> 

*/■»  %  -  (i  ^  («„-<))" 


(20) 


vd, 


a  \:w  a* 


u.  n 

tOo'l  l  -JL- ^ L - ^~L 


P/9 


q  -  q 

p  q 


Since  the  temperature  distribution  is  proportional  partly  to  *(t)  and 
partly  to  the  time  derivative  *(t),  the  perturbation  terms  H.kjm  have  two  parts 


M  •  a  #(t)  +  b..  ,  4>(t) 

lkjm  ikjm  ikjm 


(21) 
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Therefore ,  from  relation  (21)  and  the  perturbation  relation  (20),  the  re  - 
t.itivr  frequency  shift  becomes  after  integration  over  the  active  surface  of  tie 


plate, 


— -  -  a$(t)  +  ai>(t)  (22) 

o 


where  the  a^-term  is  related  to  the  first  order  temperature  coefficient  of  the 
frequency  (static  thermal  behavior)  and  a  is  the  dynamic  temperature  coefficiei: 
Theoretical  results  are  presented  in  table  I  for  different  quartz-cuts. 


Table  I 

Values  of  the  dynamic  thermal  coefficient  a  for  different  orientations 

of  quartz-crystal 


Crystal  Cuts 

ST.X 

AT  ,X 

Y.X 

((>  -  6°20' 

FST  0  -  -41“30' 

-  26° 

a  in  10~6  S/°K 

0.04 

-1.4 

-3.5 

2.8 

1 

-35.0 

-87.5 

•>*4 

o 

o 

1 

TEMPERATURE  CONTROLLED  OVEN 


SAW  delay  lines  built  with  different  cuts  of  quartz  crystal  were  submit¬ 
ted  to  temperature  cyclings  by  using  a  temperature  controlled  oven  .  As  show 
in  figure  2,  SAW  devices  and  a  LC-cut  quartz  temperature  probe  are  mounted  in 
the  same  oven.  Temperature  is  measured  with  a  quartz  thermometer  and  controlled 
by  a  computer  which  drives  the  oven  by  comparing  the  measured  temperature  with 
a  reference  temperature.  The  temperature  follows  a  variation  law  given  by  a  pn 
gram  included  in  the  computer.  Sinusoidal  cyclings  can  be  realized  with  an  acc  ■ 
racy  of  a  few  0.001°C,  as  the  SAW  oscillator  frequency  is  measured. 


Figure  2.  Temperature  controlled  oven. 

D:  SAW  delay  line;  H:  Heater; 
TP:  Temperature  probe; 

A:  Amplifier;  C:  Directional 
coupler;  DC;  D.C.  power  supply; 
TH:  Thermometer;  FC:  Frequency 
counter;  COMP:  Computer; 

PR:  Printer;  PL:  Plotter; 

FR  Freezer. 


KXPKHIMKfJTAL  BKS'dl.TS 

Figure  3  shows  an  experimental  curve  obtained  by  sweeping  temperature  o v 
a  temperature  range  of  4°C  about  15°C  with  sweep  frequency  of  2.3.10  Hz  for 

a  reguiar  ST-cut  fixed  in  two  points  by  epoxy  on  a  flat  heater.  The  measured 

_5 

dynamic  temperature  coefficient  is  equal  to  +2.10  S/°K  and  is  significantly 

larger  than  the  theoretical  one.  The  oscillator  frequency  was  87.8  MHz.  Ine  sar 

plate  was  submitted  to  temperature  cyclings  applied  on  one  end  by  a  second  heal 

source.  Thus  more  inhomogeneous  temperature  distributions  are  introduced  and 

the  measured  dynamic  temperature  coefficient  obtained  is  ten  times  greater  thar 

-4 

the  previous  one  (-2.10  S/°K). 


13  14  15  16  17 

Temperature  (®C) 


Fig.  3  -Experimental  sinusoidal  temperature  cycling.  Sweep 
frequency:  2.310-<4  Hz,  ST-Cut  AO:  +  2®K  about 
15°C  -  f13„c  -  87.8  MHz. 

- Static  frequency-temperature  characteristic. 
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CONCLUSION 


The  value  of  the  dynamic  T.C.  a,  obtained  from  the  experimental  frequen- 
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cy  responr.es  is  equal  to  2.10  s/°K  for  ST-cut  quartz  SAW  oscillator.  This  cof 
ficient  is  comparable  to  the  "a  of  AT-cut  quartz  BAW  oscillator  (1.3.10  s/°K) 

But  this  dynamic  thermal  behavior  can  be  modified  by  changing  the  mounting  of 
the  substrate  in  order  to  minimize  the  stress  effects. 

The  theoretical  model  must  be  improved  : 

-  by  introducing  the  heat  exchange  by  radiation 

-  by  using  a  bidimensional  model  which  will  take  into  account  the  edge 
cond i t ions . 

This  first  study  enabled  to  determine  the  order  of  magnitude  of  the  SAW 
dynamic  thermal  effect  and  to  make  comparison  with  bulk  waves. 

The  continuation  will  consist  in  improving  the  theoretical  model  to  study 
the  influence  of  crystal  orientation  and  wave  propagation  direction  on  the  a 
coefficient  and  determine  if  doubly  rotated  cuts  can  be  used  with  minimized  dy 
namic  effect. 
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